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We study the evolution of social networks that contain both friendly and unfriendly pairwise links 
between individual nodes. The network is endowed with dynamics in which the sense of a link in 
an imbalanced triad — a triangular loop with 1 or 3 unfriendly links — is reversed to make the triad 
balanced. With this dynamics, an infinite network undergoes a dynamic phase transition from a 
steady state to "paradise" — all links are friendly — as the propensity p for friendly links in an update 
event passes through 1/2. A finite network always falls into a socially-balanced absorbing state 
where no imbalanced triads remain. If the additional constraint that the number of imbalanced 
triads in the network does not increase in an update is imposed, then the network quickly reaches 
a balanced final state. 
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I. INTRODUCTION 

In this work, we investigate the role of friends and en- 
emies on the evolution of social networks. We represent 
individuals as nodes in a graph and a relationship be- 
tween individuals as a link that joins the corresponding 
nodes. To quantify a relationship, we assign the binary 
variable = ±1 to link ij, with Sij = 1 if nodes i and j 
are friends, and sij = — 1 if i and j are enemies (Fig.^l. 
A basic characterization of relationships between mutual 
acquaintances is the notion of social balance 0, Q ■ The 
triad ijk is defined as balanced if the sign of the product 
of the links in the triad SijSjkSki equals 1, while the triad 
is imbalanced otherwise. We define a triad to be of type 
A/j if it contains fc negative links. Thus Ao and A2 are 
balanced, while Ai and A3 are imbalanced. A balanced 
triad therefore fulfills the adage 

• a friend of my friend is my friend; 

• an enemy of my friend is my enemy; 

• a friend of my enemy is my enemy; 

• an enemy of my enemy is my friend, 

while an imbalanced triad is analogous to a frustrated 
plaquette in a random magnet Q. 

A network is balanced if each constituent triad is bal- 
anced 0, Q . An ostensibly more general definition of a 
balanced network is that every cycle in the network is bal- 
anced. Cartwright and Harary showed £| that the cycle- 
based and triad-based definitions of balance are equiv- 
alent on complete graphs. Their result implies that if 
an imbalanced cycle of any length exists in a complete 
graph, an imbalanced triad also exists. 

Balance theory has been initiated by Heider [l| and 
other social psychologists |H,|c|, and the subject remains 
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an active research area SH EHH E3- Much 
of this work was devoted to classifying balanced stable 
states of networks when relationships are static. A funda- 
mental result from these studies is that balanced societies 
are remarkably simple: either all individuals are mutual 
friends (we call such a state "paradise" ) , or the network 
segregates into two antagonistic cliques where individu- 
als within the same clique are mutual friends and indi- 
viduals from distinct cliques are enemies (we call such 
a state "bipolar") 4]. Balance theory also has natural 
applications to international relations As a par- 

ticularly compelling example, the Triple Alliance (1882) 
pitted Germany, Austria-Hungary, and Italy against the 
Triple Entente (1907) countries of Britain, France, and 
Russia 01- This bipolar state of competing alliances 
clearly contributed to the onset of World War I. 

A large network is almost surely imbalanced if the re- 
lationships are randomly chosen to be friendly or un- 
friendly. Clearly such a network is socially unstable and 
the web of relations must evolve to a more stable state 
if the individual nodes behave rationally. In this work, 
we go beyond a static description social relations and 
investigate how an initially imbalanced society achieves 
balance by endowing a network with a prototypical so- 
cial dynamics that reflects the natural human tendency 
to reduce imbalanced triads. A related line of investiga- 
tion, based on evolving social networks with continuous 
interaction strengths, has also recently appeared |15|. 



II. MODELS 

We first consider what we term local triad dynamics 
(LTD). In an update step of LTD, we first choose a triad 
at random. If this triad is balanced (Ao or A2), no evo- 
lution occurs. If the triad is imbalanced (Ai or A3), we 
change s — > — s on one of the links as follows: Ai — * Ao 
occurs with probability p, and Ai — * A2 occurs with 
probability 1 — p, while A3 — > A2 occurs with probabil- 
ity 1 (Fig. . One unit of time is defined as L update 
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events, where L is the total number of links. Notice that 
for the special case of p = 1/3, each link of an imbalanced 
triad is flipped equiprobably. When p > 1/3, the density 
of friendly links consequently tends to increase and the 
society is predisposed to tranquility, while for p < 1/3 
the societal predisposition is hostility. 

After an update step in LTD, the imbalanced target 
triad becomes balanced, but other balanced triads that 
share a link with this target may become imbalanced. 
These triads can subsequently evolve and return to bal- 
ance, leading to new imbalanced triads. Such an interac- 
tion cascade is familiar in social settings. For example, 
if a married couple breaks up, the acquaintances of the 
couple may then be obliged to redefine their relations 
with each partner in the couple so as to maintain bal- 
anced triads. These redefinitions, may lead to additional 
relationship shifts, etc. 




FIG. 1: Imbalanced triads Aj (left) and A3 (right) and the 
possible outcomes after an update step by local triad dynam- 
ics. Full and dashed lines represent friendly (e.g., sy = 1) 



and unfriendly (e. 



9; Sik 



-1) relations respectively. 



For p < 1/2, we shall show that LTD quickly drives 
an infinite network to a quasi-stationary dynamic state 
where global characteristics, such as the densities of 
friendly relations or imbalanced triads, fluctuate around 
stationary values. As p passes through a critical value of 
1/2, the network undergoes a phase transition to a par- 
adise state where no unfriendly relations remain. On the 
other hand, a finite network always reaches a balanced 
state. For p < 1/2, this balanced state is bipolar and the 
time to reach this state scales faster than exponentially 
with network size. Forp > 1/2, the final state is paradise. 
The time to reach this state scales algebraically with N 
when p — 1/2, and logarithmically in N for p > 1/2. 

We also investigate constrained triad dynamics (CTD). 
Here, we select a random link, and change s — > — s for this 
link if the total number of imbalanced triads decreases. If 
the total number of imbalanced triads is conserved in an 
update, then the update occurs with probability 1/2. Up- 
dates that would increase the total number of imbalanced 
triads are not allowed. We again define the unit of time as 
L update events, so that on average each link is changed 
once in unit of time. The global constraint accounts for 
the socially-plausible feature that an agent considers all 
of its mutual acquaintances before deciding to change 
the character of a relationship. CTD also corresponds to 
an Ising model with a three-spin interaction between the 
links of a triad, "K — — J^ijk s ij s jkSki, where the sum 
is over all triads ijk, with zero-temperature Glauber dy- 
namics As we shall see, a crucial outcome of CTD 



is that a network is quickly driven to a balanced state in 
a time that scales as In A. 

In the following two sections we analyze the dynamics 
of networks that evolve by LTD or CTD. For simplicity, 
we consider networks with a complete graph topology — 
everyone knows everyone else. This limit is appropriate 
for small networks, such as the diplomatic relations of 
countries. We then summarize and discuss some practical 
implications of our results in Sec. IV. 



III. LOCAL TRIAD DYNAMICS 

A. Evolution Equations 

We begin with essential preliminaries for writing the 
governing equations for the various triad densities. Let 
A, L — ( 2 ), and Aa = ( 3 ) be the numbers of nodes, 
links, and triads in the network. Define Nk as the number 
of triads that contain k negative links, with rik = Nk/N& 
the respective triad densities, and L + (L~) the number 
of positive (negative) links. The number of triads and 
links are related by 



_ + _3Ao + 2jV 1 +A 2 
N-2 



Ni + 2A 2 + 37V 3 
N-2 



(1) 



The numerator counts the number of positive links in all 
triads while the denominator appears because each link 
is counted N — 2 times. The density of positive links 
is therefore p = L + /L = (3uq + 2ni + n>2)/3, while the 
density of negative links isl — p = L~ /L. 

A fundamental network characteristic is N£ , which is 
defined as follows: for each positive link, count the num- 
ber of triads of type A& that are attached to this link. 
Then is the average number of such triads over all 
positive links. This number is 



(3 - k)N k 
L+ ' 



(2) 



The factor (3 — k)Nk accounts for the fact that each of 
the Nk triads of type A*, is attached to 3 — k positive 
links; dividing by L + then gives the average number of 
such triads. Analogously, we introduce N^ = kNk/L~. 
Since the total number of triads attached to any given 
link equals N — 2, the corresponding triad densities are 



K 

N-2 3n + 2ni + n 2 



(3 - k)n k 



Nr 



kn k 



N-2 m + 2n 2 + 3n 3 



(3a) 
(3b) 



We now write the rate equations that account for 
changes in the various triad densities in a single update 
event. We choose a triad at random; if it is imbalanced 
(Ai or A3) we change one of its links as shown in Fig.^ 
Let ir + be the probability that a link changes from + to 
— in an update event, and vice versa for 7r~. Since a 
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link changes from 1 — ► — 1 with probability 1 — p when 
Ai — * A 2 , while a link changes from — 1 — ► 1 with prob- 
ability p if Ai — > A and with probability 1 if A3 — ► A 2 , 
we have (see Fig. ^ 



7T + = (1 — p)ll\ 7T = pill 



n 3 . 



(4) 



Since each update changes N— 2 triads, and we also de- 
fined one time step as L update events, the rate equations 
for the triad densities have the size-independent form 



driQ 

~dt 

dni 

~dt 

dn 2 

~dT 

dn 3 

dt 



7T n n — 7T Tin 



(5) 



= 7T TIA — 7T n 
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FIG. 2: Exact stationary densities rik(p) and the density of 
friendly relations poo as a function of p. Simulation results 
for poo for N = 64 (crosses) and 256 (boxes) are also shown. 



B. Stationary States 

We first study stationary states. Setting the left-hand 
sides of Eqs. JHJ to zero and also imposing ir + = tt~ to 
ensure a fixed friendship density, we obtain 

= n>i , rii — n,2 , n,2 = . 

By forming products such as ri^n^ = n^?i7, these rela- 
tions are equivalent to 



3nori2 = rii , 3rti7i3 = n 2 



(6) 



Substituting n + and 7r _ from Eq. Q into the station- 
arity condition, tt + = ir~ , gives n 3 = (1 — 2p)ni. Using 
this result, as well as the normalization X) n fc = 1; m 
Eqs. ©, we find, after some straightforward algebra 



where 



1/[V3(1 - 2p) + 1] 



P<l/2; 
P>l/2; 



(7) 



(8) 



is the stationary density of friendly links. Equation Q 
shows that relationships are uncorrelated in the station- 
ary state. As shown in Fig. |2 the density of friendly 
links poo monotonically increases with p for < p < 1/2, 
while for p > 1/2, paradise is reached where all peo- 
ple are friends. Near the phase transition, the density 
of unfriendly links 1 — vanishes as \/3e + 0(e), as 
e = 1 - 2p 0. 



C. Temporal Evolution 

A remarkable feature of equations (0 is that if the 
initial triad densities are given by Eq. J7J — uncorrelated 



densities — the densities will remain uncorrelated forever. 
For such initial conditions it therefore suffices to study 
the time evolution of the density of friendly links p(t). 
This time evolution can be extracted from Eqs. (|3J), or it 
can be established directly by noting that p(t) increases 
if A3 — > A2 or Ai — > Ao, and decreases if Ai — > A 2 . 
Taking into account that the respective probabilities for 
these processes are 1 , p, and 1 — p, we find 



dp 
~dt 



3(2 P -i)p 2 (i-p) + (i- P y 



(9) 



Thus the time dependence of the density of friendly links 
is given by the implicit relation 



dx 



3(2p-l)x 2 (l-x) + (l-x) 3 



t. 



(10) 



When p < 1/2, the stationary density of Eq. JSJ is 
approached exponentially in time: 



Pit) 



-ct 



C 



6e 



1 



where again e = 1 — 2p. At the threshold value p = 1/2, 
the friendship density is given by 



P = l- 



1 - Po 



V'l + 2(1- A) ) 2 i 



(11) 



Here the approach to paradise is algebraic in time, viz., 
1 — p — > 1/ ' \[2~i as t — » 00. As a consequence, 



(n , ni, n 2 , n 3 ) ->• 1 = , -== , — 



1 



/2t' V2l ' 2t ' (2i) 3 / 2 
Finally when p > 1/2, 

l-p~exp[-3(2p-l)t] , 
so that paradise is reached exponentially quickly. 



(12) 



(13) 
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D. Fate of a Finite Society 

Although an infinite network reaches a dynamic steady 
state for p < 1/2, a finite network ultimately falls into 
an absorbing state for all p. Such absorbing states are 
necessarily balanced, because any network that contains 
an imbalanced triad continues to change. To see why 
such an absorbing state must eventually be reached, con- 
sider the evolution in which at each step an unfriendly 
link changes to friendly in an imbalanced triad. Since 
the number of unfriendly links always decreases, a bal- 
anced state is reached in a finite number steps. Finally, 
because this particular route to a balanced state has a 
nonzero probability to occur, any initial network ulti- 
mately reaches an absorbing balanced state. 




1 10 10 2 10 3 

N 

FIG. 3: Average absorption time as a function of N for an 
initially antagonistic society (po ~ 0) for: (a) p = 1/3; (b) 
p = 1/2; (c) p=3/4. The line in (b) has slope 4/3. 

Our simulations show that a finite network evolves to 
a bipolar state for p < 1/2, independent of the initial 
state. The size difference of the two final cliques is virtu- 
ally independent of the initial configuration because the 
network spends an enormous time in a quasi-stationary 
state before reaching the absorbing state. For p < 1/2, 
we estimate the time to reach a bipolar state by the fol- 
lowing crude argument [l7|. Consider a nearly balanced 



network. When a link is flipped on an imbalanced triad, 
then of the order of N new imbalanced triads will be cre- 
ated in the adjacent triads that contain the flipped link. 
Thus starting near a balanced state, local triad dynam- 
ics is equivalent to a biased random walk in the state 
space of all network configurations, in which the bias is 
directed away from the balanced state, with the bias ve- 
locity v proportional to N. Conversely, when the network 
is far from balance, local triad dynamics is diffusive in 
character because the number of imbalanced triads will 
change by of the order of ±iV equiprobably in a single 
update. The corresponding diffusion coefficient D is then 
proportional to N 2 . Since the total number of triads in 
a network of N nodes is ~ A 3 , we therefore expect 
that the time Tv to reach balance will scale as 

T N ~e vN */° ~e N \ (14) 

When p > 1/2, paradise is reached with a probabil- 
ity that quickly approaches one as N — > oo. At the 
threshold p = 1/2, a naive estimate for the time Tv to 
reach paradise is given by the time at which the den- 
sity of unfriendly links a(t) = 1 — p(t) is of the order 
of A~ 2 , corresponding to one unfriendly link in the net- 
work. From Eq. (|11|) . the criterion a(Tjv) ~ A~ 2 gives 
Tjv ~ A 4 . While simulations show that Tv does scales 
algebraically with A, the exponent value is much smaller 
(Fig. Bib)). The source of this smaller exponent is the 
existence of anomalously large fluctuations in the number 
of unfriendly links. 

To determine these fluctuations in the thermodynamic 
limit, we write the number of unfriendly links A(t) = 
L~(t) in the canonical form 0] 

A(t)=La(t) + VLrj{t), (15) 

where a(t) is deterministic and rj(t) is a stochastic vari- 
able. Both a and 77 are size independent in the thermo- 
dynamic limit (L ^> 1), and the form of Eq. I|15l) assures 
that the average (A) and the variance (A 2 ) — (A) 2 grow 
linearly with the total number of links L. In AppendixlAl 
we show that a = (r/ 2 ) grows as 

as t — > 00 . (16) 

Thus the time to reach paradise Tv is determined by 
the criterion that fluctuations in A become of the same 
order as the average, viz., 

v/La(Tjv) ~ La(T N ) . (17) 

Using a(t) ~ from Eq. Jill) together with Eq. JEJ 

and L - N 2 , we rewrite Eq. (T^J as N 2 T N 1/2 - iVT^/ 4 . 
This leads to the estimate 

Tv - N 4/3 . (18) 

Above the threshold p > 1/2, paradise is approached 
exponentially quickly [see Eq. (|13fl ] and the time to par- 
adise scales logarithmically with network size: 

Tv - (2p- I)- 1 In AT (19) 
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Interestingly, the estimates (|18f) and l|19f) coincide when 
2p - 1 ~ 7V-4/3 i n TV. That is, there is a finite-width 
critical region near the phase transition due to finite- 
size effects. In Appendix ^ we estimate this width by 
analyzing the fluctuations below the threshold p < 1/2 
and obtain essentially this same result. 

Summarizing, the asymptotics for the absorption time 
are: 



'exp(A 2 ) P<l/2 
T N ot{N 4 / 3 P=l/2 
(2p~l)- 1 \nN p>l/2 

in agreement with the simulation results in Fig. [3| 
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IV. CONSTRAINED TRIAD DYNAMICS 

A. Jamming and Absorption Time 

In constrained triad dynamics (CTD), the number of 
imbalanced triads cannot increase in an update event, 
and the final state can either be balanced or jammed. 
A jammed state is one in which imbalanced triads exist 
and for which the flip of any link increases the number 
of imbalanced triads. Since this type of update is for- 
bidden in CTD, there is no escape from a jammed state. 
Moreover, jammed states turn out to be much more nu- 
merous than balanced states 



(see Sec. IIV Dl . In spite 
of this fact, we find that the probability for the network 
to reach a jammed state, Pj am (A), quickly goes to zero 
as N increases, except for the case of an initially antag- 
onistic society (po — 0), where Pj am (iV) decays slowly 
with N (Fig. 0} . Thus the final network state is almost 
always balanced for large N, and consists either of one 
clique (paradise) or two antagonistic cliques. It is worth 
mentioning that we never observed "blinkers" m sim- 
ulations, although we cannot prove that such states do 
not exist. These are trajectories in the state space that 
evolve forever and correspond to a network in which all 
possible updates involve no change in the number of im- 
balanced triads. 

Another fundamental feature of CTD is that the time 
Tjv for a network of N nodes to reach its final state gener- 
ically scales as In N. While Tm now depends on the initial 
condition, in contrast to LTD, this dependence occurs ei- 
ther in the amplitude or in lower-order corrections of the 
absorption time. Thus the logarithmic growth of Tm with 
TV is a robust feature of CTD. 



B. Final Clique Sizes 

An unexpected feature of CTD is the phase transi- 
tion for the difference in sizes C\ and Ci of the two fi- 
nal cliques as a function of p (Fig. [5J . We quantify this 
asymmetry by the scaled size difference 6 = (Ci — C^/N, 
F° r Po ^ 0-4 the cliques sizes in the final bipolar state 



FIG. 4: Probability of reaching a jammed state Pj a 
function of N for several values of po- 



are nearly the same size and (S 2 ) ~ 0. As po increases 
toward p$ w 0.65, the size difference of the two cliques 
continuously increases. A sudden change occurs at pj, 
beyond which the final network state is paradise. The 
probability distribution for S is sharply peaked about its 
average value as N — > oo (Fig.[SJ). Since (S 2 ) and the den- 
sity of friendly links p^ are related by {S 2 ) = Ip^ — 1 
in a large balanced society, uncorrelated initial relations 
generically lead to p,^ > po. Thus CTD tends to drive a 
network into a friendlier final state. 



CO 
V 



0.6 



0.4 ■ 



0.2 



N = 256 o 
512 

1024 v 

2048 □ 



k <$? A Gp A Oft I 



0.2 



0.4 0.6 

Po 



FIG. 5: Asymmetry of the final state as a function of the 
initial friendship density po for several network sizes. 

While we do not have a detailed understanding of this 
phase transition, we give a qualitative argument that sug- 
gests that a large network undergoes a sudden change 
from poo = (two equal size cliques) when po < 1/2 to 
Poo = 1 (paradise) when p > 1/2. The fact that the 
transition appears to be located near p$ w 0.65 (Fig. |5j) 
rather than at po = 1/2 indicates that our approach is 
not a complete description for the transition. 

We first assume, as observed in simulations of large 
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FIG. 6: Scaled probability distribution of the relative differ- 
ence between the final clique sizes for po — and 1/2. 



networks, that jammed states do not arise. We also 
assume that a network remains uncorrelated during its 
early stages of evolution. Consequently the densities 
n + = (no , nf, n\, n%) of triads that are attached to a 



positive link are 



(p 2 ,2p(l-p),(l-p) 2 ,0). 



(21) 



to — , it is necessary that 



For a link to change from 

+ > Hq + n~2 ■ From Eq. 121fl . this condition is 
equivalent to 4p(l — p) > 1, which never holds. Similarly, 
the densities n~ = (rig , , , nj) of triads attached to 
a negative link are 



n- = (0,p 2 ,2p(l~p),(l-p) 2 ). 



(22) 



The requirement n\ + > Uq + now reduces to 
1 > 4p(l — p), which is valid when p ^ 1/2. 

Thus for a large uncorrelated network, only negative 
links flip in CTD. Since, the density of negative links is 
1 — p, the governing rate equation is 



dp 



1-P, 



from which 



p = 1 - (1 - p )e~ 



(23) 



(24) 



From the criterion 1 — p{Tn) ~ ^V~ 2 , corresponding to 
one unfriendly link remaining in the network, the time 
to reach paradise is given by Tn ~ In N , in agreement 
with simulations. According to Eq. (|24|) a network should 
evolve to paradise for any initial condition. 

However, our simulations indicate that this homoge- 
neous solution is unstable for po < 1/2. In this case, the 
density of friendly links p initially still increases according 
to l|24l) until p w 1/2. At this point, correlations in the 
relationship structure begin to arise and these ultimately 
lead to a bipolar society with p x ps 1/2. We now give a 
qualitative argument to support these observations. 



When p(i) sa 1/2, there are many partitions of the 
network into two subnetworks §i and §2 of nearly equal 
sizes C\ = I Si I and C2 = | S 2 1 , for which the density of 
friendly links within each subnetwork, p\ and P2, slightly 
exceed 1/2, while the density (3 of friendly links between 
subnetworks is slightly less than 1/2. Our basic point is 
that this small fluctuation is amplified by CTD so that 
the final state is two nearly equal-size cliques. 

To appreciate how such an evolution can occur, we 
assume that relationships within each subnetwork and 
between subnetworks are homogeneous. Consider a neg- 
ative link in Si. The densities of triads attached to this 
link are given by (|22|) . with p replaced by (3 when the 
third vertex in the triad belongs to §2 , and by l|22l) , with 
p replaced by pi when the third vertex belongs to Si. 
The requirement that a link can change from — to + 
according to CTD now becomes 

Ci[l - 4pi(l - Pl )] +C 2 [1- 4/5(1 - /?)] > 0, (25) 

which is always satisfied. Additionally, negative links 
within each subnetwork can change to positive with rate 
one, while positive links within each subnetwork can 
never change. 

Consider now a positive link between the subnetworks. 
The triad densities attached to this link are given by 

n+ = (/3p„/3(l - Pi ) + Pj (l ~ 13), (1 - - Pj ), 0) 

when the third vertex belongs to Sj . Since 

/3(l-p J )+p,(l-/3)-/3p,-(l-/3)(l-p J ) = (2p 3 -l) (1-2/3). 

the change H > — is possible if 

[Ci(2pi-1) + C 2 (2p2- l)](l-2/3) >0. (26) 

Thus if the situation arises where pi > 1/2, p 2 > 1/2, and 
f3 < 1/2, the network subsequently evolves to increase the 
density of intra-subnetwork friendly links and decrease 
the density of inter-subnetwork friendly links. These link 
densities thus evolve according to the rate equations 



dpi , dp 2 



and give the instability needed to drive the network to a 
final bipolar state. 

The last step in our argument is to note that when 
C\ ~ C2 ~ N/2, the number of ways, ( c ), to partition 
the original network into the two nascent subnetworks Si 
and S2, is maximal. Consequently, the partition C\ = C2 
has the highest likelihood of providing the initial link den- 
sity fluctuation, after which the homogeneous evolution 
(|23|l is replaced by the clique evolution l|27l) so that a 
homogeneous network organizes into two nearly equal- 
size cliques. Although our argument fails to account for 
the quantitative details of the transition shown in Fig. [5] 
the primary behaviors of (S 2 ) in the two limiting cases of 
Po — > and po — > 1 are described correctly. 
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C. Structure of Jammed Configurations 

While jammed configurations can arise in CTD, wc will 
now show that: Jammed states are possible if and only 
if the network size is TV — 9 or TV > 11. To prove this 
statement, we first explicitly construct jammed config- 
urations for N — 9 and TV > 11. Fig. [7| shows three 
jammed configurations for TV = 9, the smallest possible 
TV where jammed configurations can occur. The exam- 
ple in Fig. 13 a) was observed in simulations, while the 
jammed configuration in Fig. 0[b) consists of three an- 
tagonistic cliques of 3 nodes each. We now generalize this 
latter construction of jammed states to arbitrary TV > 11. 

Consider three mutually antagonistic cliques of sizes 
(mi, rri2,ms), with mi + m 2 + 7723 = TV. A link within 
a clique is necessarily stable, as all attached triads are 
of type Ao or A2. Conversely, a negative link between 
clique 1 (circles in Fig. [3b)) and clique 2 (squares) is 
attached to both stable and imbalanced triads. There 
are mi — 1 + 7772 — 1 attached stable triads of type A2, 
where the third node of the triad is either within clique 
1 or clique 2, and 7713 attached imbalanced triads of type 
A3, where the third node is in clique 3 (triangles). The 
requirement for link stability among these cliques is then 

mi + 7772 > m,3 + 2 

m 2 + 7773 > mi + 2 (28) 

7773 + mi > m 2 + 2, 

where the last two equations arise by cyclic permutations 
of the first. 

We term a partition (mi, 777,2, 777,3) "jammed" if it satis- 
fies the inequalities l|28|l . By summing pairs of Eqs. Q28JI. 
we find 777j >3,j= 1, 2, 3. Thus jammed partitions are 
possible only for networks of size TV = mi + 7772 +777,3 > 9. 
Following the rules in Eq. 1281) . the following partitions 
are jammed: for N — 3k with k > 3, partitions of the 
form (k,k,k); for TV = 3fc + 2 with k > 3, (fc,fc+l,fc+l); 
and for TV = 3k+l with k > 4, (fc, k, fc+1). Thus jammed 
partitions indeed exist for TV = 9 and TV > 11. 

Finally, we show that jammed states are impossible 
for TV < 8 and TV = 10. As a preliminary, we need the 
following: 

Lemma. Let ABC be imbalanced, ADC balanced, and 
sac = — 1. Then one of the two triads ABD and BDC 
is balanced, and the other is imbalanced. 

Proof. Let sxyz — sxySyzszx be the sign of the 
triad XYZ. For the imbalanced ABC triad we have 
sabc = ~~ 1 while for the balanced triad sacd = L Using 
additionally the identities sxy — syx, s XY = 1, we 
obtain the product of the signs of triads ABD and BDC: 

sabdSbdc = sabsbdsdasbdsdcscb 

= sabsdasdcscb 

= sabsbcscasacscdsda 

— sabcsacd = -1 

from which the lemma immediately follows. 




(b) (c) 

FIG. 7: Examples of jammed configurations for TV = 9 (only 
friendly links are displayed), (a) A jammed configuration that 
appeared in simulations, (b) A jammed state consisting of 
three mutually antagonistic cliques, (c) A jammed state de- 
rived from (b) in which the top clique from (b) is friendly 
toward the remaining two cliques. 



Now suppose that there is a jammed state in a net- 
work with an even number of nodes TV = 2k. By def- 
inition, there is at least one imbalanced triad in the 
jammed state; let ABC be such an imbalanced triad 
with sac = ~ 1- Since the state is stable, out of the 
TV — 2 = 2k — 2 triads attached to the link AC, at least 
k are balanced. (Note that this construction requires 
k > 3; however, it is trivial to show that there is no 
jammed state for the case k — 2 [TV = 4]). Take k such 
balanced triads and denote them ADjC, j — 1, . . . , k. 
To each pair of triads ABC and ADjC we now apply the 
lemma. Then there is a certain number x of imbalanced 
triads among ABDj, and a certain number y of imbal- 
anced triads among CBDj, with x + y = k. Stability 
ensures that there are at most k — 2 imbalanced triads 
attached to the link AB. Recalling that ABC is imbal- 
anced and that there are x imbalanced triads ABDj , we 
obtain x + 1 < k — 2. A similar argument applied to link 
CB leads to y + 1 < k — 2. Summing these inequalities 
and using x + y = k gives k > 6, or N > 12 for even TV. 

The case of odd TV is similar. We set TV = 2k — 1, with 
k > 3. Now each link is attached to at least k balanced 
triads and at most k — 1 imbalanced triads. Repeating 
the same argument as for the even case we obtain the 
conditions x + 1 < k — 1 and y + 1 < k — 1 with x-\-y = k, 
which leads to N > 9 for odd TV. 
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D. Number of Jammed Configurations 



SUMMARY AND DISCUSSION 



Last, we show that the number of jammed configu- 
rations greatly exceeds the number of balanced config- 
urations. The total number of distinct network con- 
figurations is 2 L . Each balanced state has the form 
{m\,m2\m\ + m 2 = N), and we enumerate all classes 
of balanced states by counting the integer solutions of 
mi + m 2 = N , with < m% < m%. Therefore the num- 
ber of classes of balanced states is B = k + 1 for N = 2k 
and N = 2k + 1. The total number 23 of balanced states 
is determined from 



23 



E 

mi+m2- 



mi + m 2 
m\ 



(29) 



and is thus much larger than the number of classes of 
balanced states. 

For the number of classes of jammed states J and the 
number of jammed states 3, we can only establish lower 
bounds. For large N, instead of exact counting we em- 
ploy a continuum description. From Eq. I|28[) . the number 
of jammed partitions is equal to N 2 times the area A of 
the region inside the triangle xi + x 2 + £3 < 1 defined by 
inequalities x\ + x 2 > £3, ^2+^3 > Xi, £3 + x\ > x 2 ; 
this area is A = 1/8. We also divide by 3! = 6 to account 
for over-counting different permutations of mi, 7712,7713. 
Thus J > A 2 /48. We could improve this bound by 
counting additional jammed states built from the con- 
struction in Fig.|7fc), but this contribution would not af- 
fect the N dependence of the bound. However, we do not 
know whether J oc iV 2 or J grows faster than N 2 due to 
the existence of jammed states, such as those in Fig. Ufa), 
that are not in the classes described in Sec. IIV CI 

We obtain a lower bound for 3 in a similar manner 
to that in Eq. I)29|l for counting the number of balanced 
states, 



3> E 

jammed 



mi+ m 2 + m 3 
mi, m 2 , 777.3 



(30) 



where the sum is over jammed partitions, and the trino- 
mial coefficient is 



777l + 777 2 + 7773 
777l, 777 2 , 777 3 



(777,1 + 777 2 + 77l 3 )! 
7771 ! 7772! 7773! 



The summand in Eq. (|30[l is sharply peaked around mi = 
7772 = 7773 = JV/3 and therefore the sum is very close 
to 3 N which is the sum over all partitions. Again, the 
lower bound may be weak because of the neglect of non- 
tripartite jammed configurations. 

In summary, we find that 3 > 3 N > 2^ = 23. Thus the 
total number of jammed states greatly exceeds the total 
number of balanced states. Nevertheless, for a random 
initial condition, the probability to end in a balanced 
state is very close to one while the probability to end in a 
jammed state is negligible; that is, the basin of attraction 
of balanced states greatly exceeds the basin of attraction 
of jammed states. 



In social relations, we may encounter the uncomfort- 
able situation of an imbalanced triad. If you have two 
friends that develop a mutual animosity, then an imbal- 
anced triad of relations exists. You will then likely have 
to choose between these two friends, thereby resolving 
the social conflict and restoring the relationship triads to 
balance. In this work, we implemented simple and pro- 
totypical dynamical rules for healing imbalanced triads 
and we investigated the resulting evolution of these social 
networks. 

In the case of local triad dynamics, a finite network 
falls into a socially-balanced state, where no frustrated 
triads remain. The time to reach this final state depends 
very sensitively on the propensity p for forming friendly 
links in the update events that heal social imbalance. For 
an infinite network, the balanced state is never reached 
when p < 1/2 and the system remains in a stationary 
state. The density of unfriendly links gradually decreases 
and the network undergoes a dynamical phase transition 
to an absorbing, paradise state for p > 1/2. 

We also examined the dynamics in which an additional 
global constraint is imposed that the number of imbal- 
anced triads in the entire network cannot increase in an 
update event. The virtue of this dynamics is that the 
final outcome is always reached quickly. A downside, 
however, is that the final configuration of the network 
may be jammed — these are states that are not balanced, 
but where flipping any link increases the number of im- 
balanced triads. Fortunately, the probability of reaching 
a jammed state is vanishingly small and the final state is 
either a two-clique bipolar state or paradise. 

As alluded to in the introduction, a natural applica- 
tion for social balance ideas is to international relations, 
with the prelude to World War I being a particularly 
appropriate example. For example, the Three Emper- 
ors' League (1872, and revived in 1881) aligned Ger- 
many, Austria-Hungary, Russia, leaving France isolated. 
However, the Turkish-Russian war (1877) and tension 
between Austria-Hungary and the Balkan states unrav- 
eled Russia's participation in the League, and a bipartite 
agreement between Germany and Russia lapsed in 1890. 
In the meantime, the Triple Alliance was formed in 1882 
that joined Germany, Austria-Hungary, and Italy into a 
bloc that continued until World War I. 

On the other hand, a French-Russian alliance was 
formed over the period 1891-94 that ended France's diplo- 
matic isolation with respect to the Triple Alliance. Sub- 
sequently an Entente Cordiale between France and Great 
Britain was consummated in 1904, and then a British- 
Russian agreement in 1907, after long-standing tensions 
between these two countries, that then bound France, 
Great Britain, and Russia into the Triple Entente. While 
our historical account of these Byzantine maneuvers is 
very incomplete (see Refs. |l4j for more information), the 
basic point is that among the six countries that comprised 
the two major alliances, bipartite relationships changed 
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as triads became unbalanced and there was a reorganiza- 
tion into a balanced state of the Triple Alliance and the 
Triple Entente that became the two main protagonists at 
the start of World War I. 

On the theoretical side, there are several avenues for 
additional research. One possibility is to relax the defi- 
nition of imbalanced somewhat. This is the direction fol- 
lowed by Davis [|j who proposed the "clusters model" in 
which triads with three unfriendly relations are deemed 
acceptable. The clusters model thus allows for the pos- 
sibility that "an enemy of my enemy may still be my 
enemy." This more relaxed definition for imbalanced tri- 
ads may lead to interesting dynamical behavior that will 
be worthwhile to explore. 

Another natural generalization of the balance model 
would be to ternary relationships of positive +, nega- 
tive — , or indifferent 0. These relations may lead to the 
emergence of cliques (groups of mutual friends who dis- 
like other people) and communities (groups of mutual 
friends with no relations with other people). It would be 
interesting to study the number of cliques and number of 
communities as a function of network size and the density 
of indifferent relationships. Communities on the Web can 
be effectively identified 20] , and these results may allow 
for useful comparisons between data and model predic- 
tions. Finally, relations need not be symmetric, that is, 
Sij may be different from Sji , and it may be interesting to 
generalize the basic notions of balance to networks with 
such asymmetric interactions. 



Acknowledgments 

TA gratefully acknowledges financial support from the 
Swiss National Science Foundation under the fellowship 
8220-067591. SR thanks K. Kulakowski for informative 
discussions about social balance, and also acknowledges 
financial support from DOE grant W-7405-ENG-36 (at 
LANL) and NSF grant DMR0227670 (at BU). 



Ai — > N2, respectively. From (|A1|) we obtain 



-{A) = -{N 3 )-p{N 1 ) + (l-p){N 1 }, (A2) 



which simplifies to 



dt 



(A) = -(A 3 ) 



(A3) 



at the threshold p = 1/2. Since (A) cx a and (A3) cx 
(A 3 ) cx a 3 to lowest order, Eq. (|A3|) can be written as 



da 
~dt 



(A4) 



whose solution is given by Eq. Qlip. Similarly from (|A1|) 
we obtain 

at 

+ {1- P )({2A+1)N 1 ) (A5) 
which, for p = 1/2, simplifies to 
d 



dt 



(A 2 ) = (A 3 ) + (Ai) - 2 (A A 3 ) 



(A6) 



Subtracting Eq. (|A3|) multiplied by 2(A) from Eq. JXBJ, 
we obtain the evolution equation for the variance 

I ((A 2 ) (A) 2 ) = (A 3 ) + (Ax) + 2 ((A)(N 3 ) - (A A 3 )) . 

Now using standard methods 0] to compute moments 
of the stochastic variable A = La + \TLr\ <|15fl , we obtain 
the leading behavior (A 2 ) — (A) 2 cx a, while (A3) cx a 3 , 
and (Ai) cx 3a (1 — a) 2 . Similarly, the leading terms in 
(A) (A3) and (A A3) cancel, while the next correction is 

(A) (A 3 ) - (A A3) cx 3a 2 CT - 6a 2 o- . (A7) 

Using these results for the various moments, the variance 
satisfies 



APPENDIX A: FLUCTUATIONS IN LOCAL 
TRIAD DYNAMICS 

In this appendix, we compute the normalized variance 
a = (r/ 2 ). Wc focus on the most interesting case of the 
critical regime p = 1/2, where fluctuations exhibit the 
asymptotic behavior of Eq. l|16fl . Then we briefly discuss 
the two regimes p < 1/2 and p > 1/2. 

Wc first note that A changes according to 



A - 1 rate A 3 

A - 1 rate pNi 

A + l rate (l-p)Ai 



(Al) 



which describe the processes A3 — > A2, Ai — > Aq, and 



— = -6a 2 ( 7 + a 3 + 3a(l- a ) 2 
dt 

Dividing (|A8|) by (|A4|) we obtain 



da _ 6 a 3 + 3a(l-a) 2 
da a a 3 



(A8) 



(A9) 



Since a = a 6 solves the homogeneous equation da /da = 
6a /a, we seek a solution of the inhomogeneous equation 
i9\\ in the form a = a 6 s(a). Equation (|A9I) becomes 



ds 
da 



a 3 + 3a(l - a) 2 



3 6 4 

a 8 a 7 a 6 



whose solution is 



3 1 1 4 1 ^ 

7 a 7 a 6 5 a 5 



(A10) 
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Thus a = a 6 s(a), with s(a) given by (|A10() . The inte- 
gration constant C is fixed to satisfy the initial condition 
<r(ao) = 0. In particular, for a totally antagonistic initial 
network (p — or a = 1), C = —8/35. In this case 
a 



1 / y/1 + 2t, and the variance becomes 



3 1 
7 a 



35 



(All) 



The leading asymptotic behavior a — > (i\f2/7)t 1 / 2 
holds independent of the initial condition. Hence we es- 
tablish the crucial result l|16l) , which leads to the asymp- 
totic behavior l|18|) for the absorption time for p = 1/2. 

For p ^ 1/2, or e = 1 - 2p 0, we recast Eq. (jA2() to 



da , s0 

— = 3eo(l - a? 



(A12) 



which is of course identical to Eq. ©. Following the 
same steps that led to Eq. l|A8ll . we then derive for the 
variance 

^ = a 3 + 3a(l - a) 2 - 6a[a 2 - e(l - a)(l - 3a)] . (A13) 

When p < 1/2, both a and er quickly approach station- 
ary values 



'lie 



/3e+ 1 



1 



1 



(A14) 



The density of unfriendly links in a finite system there- 
fore exhibits fluctuations of the order of \J o^jL about 
the average density a^. Close to the phase transition 
point, the magnitude of fluctuations eventually becomes 
comparable with the average. From ~ U o^jL 
and Eq. (| A14|) . we find that this equality occurs when 
e ~ A^~ 4 / 3 ; this gives an estimate of the width of the 
phase transition region due to finite-size effects. 

When p > 1/2, both a and a vanish as t — ► oo. A 
straightforward asymptotic analysis of Eq. I|A13|I yields 



a 



2p- 1 



as a 



(A15) 



Fluctuations become comparable with the deterministic 
part when a ~ a jL ~ */ a/(2p — 1)L, that is, when 
La ~ {2p - l)" 1 . Using a - e ~^v-^ f rom Eq. (O we 
estimate the time to reach paradise to be 



T JV ~(2p-l)- 1 ln[(2p-l)Af 2 ]. 



(A16) 



The difference between this result and Eq. (|19fl . which 
was established using the naive criterion La ~ 1, is small 
because the factor 2p — 1 appears inside the logarithm. 
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